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Abstract. Density functional theory (DFT) calculations of the energy of magnetic 
anisotropy for diluted ferromagnetic semiconductor Gei_a,Mn 2 ,Te were performed 
using OpenMX package with fully relativistic pseudopotentials. The influence of 
hole concentration and magnetic ion neighbourhood on magnetic anisotropy energy 
is presented. Analysis of microscopic mechanism of magnetic anisotropy is provided, 
in particular the role of spin-orbit coupling, spin polarization and spatial changes 
of electron density are discussed. The calculations are in accordance with the 
experimental observation of perpendicular magnetic anisotropy in rhombohedral 
Gei_a;Mn 3 ,Te (111) thin layers. 
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1. Introduction 

Microscopic understanding of anisotropic properties of any ferromagnet constitutes a 
very challenging theoretical task and is very important from the point of view of 
possible applications. Particularly, in the case of ferromagnetic semiconductors, a deep 
understanding of the mechanisms responsible for the behavior of magnetization becomes 
crucial because a number of semiconductor properties may be changed by external 
factors. It was established that direction of the easy axis of magnetization can be 
changed by changing the density of free carriers or by applying a strain. The carrier 
concentration can be varied by applying gate voltage to the layer, and strain in a him 
can be imposed by cementing it to a piezoelectric actuator. These and other related 
phenomena are described in a recent review pp. 

Previous theoretical descriptions of magnetic anisotropy in semiconductors were 
based either on the mean-held Zener model, or on the concept of single ion anisotropy. 
The Zener model was applied to III-V and II-VI diluted ferromagnetic semiconductors 
[21E]. In this approach, it is assumed that the nonzero macroscopic magnetization of 
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magnetic ions leads, via sp — d exchange interaction, to the spin splitting of band states. 
Because of the presence of spin-orbit (S-0) interaction, the band structure, and thus the 
energy of carriers populating band states, depend on the direction of magnetization. In 
this model the magnetic anisotropy is due to the simultaneous presence of two factors, 
i.e. the hnite concentration of free carriers, and the non-vanishing S-0 coupling in 
the host crystal. For insulating crystals, or for the vanishing band S-0 interaction, a 
vanishing magnetic anisotropy is predicted. (Note that the limit of isolated impurities is 
beyond the scope of the model.) The Zener model explained, at least semiquant it avely, 
a number of experiments [H U El E]. 

The second approach is based on the concept of single ion anisotropy. The anisotropic 
properties of a macroscopic crystal are brought about by anisotropic magnetic properties 
of magnetic 3d or 4/ ions subject to ligand electric helds of the crystal environment. 
In II-VI and IV-VI semiconductors, the Mn^’*' ions in the crystal retain to some extend 
their atomic-like character. In the absence of the S-0 coupling, according to the Hund’s 
rule hve electrons on 3d shell form an orbital singlet with the total angular momentum 
L = 0. This system does not interact with the crystal held, and the ground state 
of the ion is sixfold degenerate since the total spin S = b/2. Actually, however, the 
Hund’s rule is an approximation only, because due to relativistic S-0 effects there is 
a coupling between and the excited states of 3d^ conhguration. Consequently, 
the ground state is no longer a pure S state, and the coupling of Mn with its crystal 
surrounding becomes possible. This leads to the ground state splitting, and Mn can be 
described by the 6x6 matrix of the effective spin hamiltonian [7]. For example, the spin 
hamiltonian of Mn^"*" in a cubic environment is proportional to Sf. + Sy + S^. This effect 
is experimentally observed as a hne structure splitting of the electron paramagnetic 
resonance (EPR) spectra [^. During the last hve decades various calculations of this 
ehect were presented, see for example Refs. ii El [iDi [n]. However, in these works the 
virtual, excited states taken into account in the perturbation treatment had the same 
number of electrons on 3d orbitals as the ground state, i.e., hve. 

This constraint was lifted in another approach to the single ion anisotropy proposed 
in Ref. [12] . In this work, the system consisted of the Mn^"*" ion and band carriers, and in 
the perturbation calculations of the ehective spin hamiltonian the number of electrons 
on the excited 3d states could be diherent than hve. In other words, the virtual transfer 
of electrons between the 3d shell and band states was allowed. This idea was previously 
widely used in the calculations of sp — d kinetic exchange integral in semiconductors 
[131 HD, or in the calculations of exchange coupling between magnetic ions [T3] . 

Within this approach it was possible to study the inhuence of both carrier 
concentration and alloy disorder on the ground state splitting. The method succesfully 
explained the EPR measurements of Mn ion in strained PbTe layers [T2|, the single 
ion magnetic anisotropy in Pbi_a;_j/Sna;MnyTe mixed crystals [I6|, the magnetic specihc 
heat in Pbi-ajMn^-Te [T7], as well as that of PbTe doped with Eu^+ and Gd^+ |18j . 

In the above theories of single ion magnetic anisotropy the manganese spin is treated 
as a quantum mechanical object. In particular, the nonvanishing S-0 coupling for 
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3d electrons is necessary to prodnce a nonzero gronnd state splitting, or single ion 
anisotropy of Mn. The band S-0 conpling is of minor importance in the sense that 
even for the zero band S-0 interaction we get the splitting, althongh much smaller 
than in the case of the actnal S-0 conpling. On the other hand, it is also possible to 
propose a model of single ion anisotropy [I9] in which the manganese spin, like in the 
Zener model, is treated as a classical Heisenberg vector interacting with band carriers 
via sp — d interaction, bnt of conrse, in this case the nonzero band S-0 conpling is 
necessary. In this model, a strong dependence of magnetic anisotropy constants on the 
concentration of free holes was predicted and qualitatively confirmed experimentally [20] 
in Sni_ 2 ,Mnj;Te. 

In the present paper we study magnetic anisotropy in Gei_ 2 ;Mna,Te ferromagnetic 
semiconductor by first principle calculations. A natural question arises what, concerning 
the magnetic anisotropy, significantly new information may be obtained applying this 
method? Certainly, using DFT calculations it is not possible to obtain directly the 
ground state splitting or the effective spin hamiltonian for the Mn ion. This is not related 
to the accuracy problem or the limitations in computational resources but to the fact 
that DFT calculations determine the charge density in the ground state of the system. 
On the other hand, as it is shown below, this approach allows to investigate microscopic 
mechanisms responsible for the anisotropy, in particular the role of S-0 coupling, of the 
range of spin polarization of atoms surrounding the magnetic ion, and of the influence 
of crystal disorder. Although we concentrate exclusively on Gei_a,Mna;Te we think that 
the general conclusions may be useful for other ferromagnetic semiconductors. 

GeTe is a nonmagnetic narrow gap IV-VI semiconductor, which has rhombohedral 
crystal symmetry with the angle characterizing the unit cell a ~ 88.35° at room 
temperature. In addition, the Ge and Te sublattices are displaced relative to each other 
along the [111] direction by a vector rao(lll), where the lattice constant Oq = 5.98 A 
and r ~ 0.03 [221123]; this displacement causes GeTe to be a unique ferroelectric. These 
features also characterize magnetic Gei_a;Mn 3 .Te mixed crystals for moderate manganese 
contents x < 0.2. 

In IV-VI semiconductors, Mn ions are incorporated as Mn^’*' with the 3d^ electron 
conhguration. They introduce local magnetic moments and remain electrically 
neutral. Gei_a;Mna;Te diluted magnetic semiconductor exhibits ferromagnetic transition 
with the Gurie temperature Tc{x) up to 190 K [251 12S]- In these mixed crystals 
ferromagnetism is induced by quasi-free carriers. The simultaneous presence of 
ferromagnetic and ferroelectric order (multiferroicity) makes this compound very 
interesting for possible applications, in particular as a magnetic phase-change material, 
or a spintronic system with electrical polarization controlled spin splittings of band 
states [27| . 

Recently, magnetization and ferromagnetic resonance (FMR) measurements were 
performed on both monocrystalline layers of Gei_a,Mn 3 ;Te grown on the BaF 2 (HI) 
substrate and on polycrystalline, layered Geo.gMno.iTe microstructures [231 ESI EEl ESI 
[30] . The experiments probed magnetic anisotropy energy (MAE) in the ferromagnetic 
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state. One of the striking results is that in monocrystalline layers with the manganese 
content of the order of 10 at. %, the easy axis of magnetization is oriented along the 
[111] crystallographic direction perpendicular to the layer, while the usual in-plane 
easy axis due to dipolar interactions (shape anisotropy) is observed in polycrystalline 
microstructures and is partially recovered in strained monocrystalline Gei_a;Mn 3 ;Te 
layers upon annealing. In the layers with the Mn content higher than 20 at. %, for 
which the X-ray diffraction studies reveal the cubic (rock-salt) structure, the usual in¬ 
plane easy axis is observed. 

Motivated by the experimental results we study microscopic mechanisms of 
magnetic anisotropy and the decisive factors that determine MAE. We mainly report 
the results obtained for rhombohedral crystals characterized by Oq = 5.98 A, a = 88" 
and T ^ 0.03, and low Mn concentrations. Most of the calculations were performed 
for 2 X 2 X 2 supercells containing 64 atoms with one Ge replaced by Mn atom, which 
corresponds to approximately x = 0.03 in Gei_ 3 ;Mn 2 ,Te mixed crystal. For comparison, 
we present also some results for the 216 atoms 3x3x3 supercell containing one Mn 
atom, which corresponds to x ~ 0.01. These Mn concentrations are much smaller 
than those in the samples studied experimentally, and this is the hrst reason why a 
quantitative comparison with experiment is not possible. The second reason is also 
important: preliminary calculations clearly indicated that MAE strongly depends on 
the microscopic disorder present in an alloy. This disorder is inevitable in mixed crystals 
due to both the different ionic radii of ions constituting the crystal, and to the random 
occupation of lattice sites. Gonsequently, in an actual crystal each Mn ion is in a different 
surrounding. In general, these surroundings have no symmetry, what leads to different 
directions of the easy axis of magnetization for different Mn ions. Below, an example 
of the dependence of MAE on the nearest neighbourhood of Mn ion is given. However, 
a detailed analysis of influence of the microscopic disorder is beyond the scope of the 
paper. That is why geometry optimization of the supercells typically was neglected. 

The paper is organized as follows. After presenting the technical details of 
calculations in Section II, in Section III we analyze the dependence of MAE on the 
concentration of free holes and on the manganese ion’s neighbourhood. Section IV 
contains analysis of various factors that determine MAE. We discuss the role of spin- 
orbit coupling, of spin polarization of free carriers in the supercell, and of the changes 
of the electron density caused by the changes of manganese spin’s direction. Section V 
summarizes the paper. A few preliminary results were presented in a short conference 
communication [3T] . 

2. Technical details of calculations 

The calculations were performed with the open-source OpenMX package [32] with 
fully relativistic pseudopotentials. The Local Density Approximation together with the 
Geperly-Alder [33] exchange-correlation functional was used. The inclusion of spin-orbit 
interaction is crucial for calculation of energy of magnetic anisotropy, which disappears 
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when the spin-orbit interaction is omitted. 

2.1. Pseudopotentials 

The pseudopotentials distributed with OpenMX are generated assuming 16 valence 
electrons for Te. Because of the very high accuracy demanded by the studied problem, 
large sizes of supercells and very high density of grids in the k-space, it was necessary 
to use pseudopotentials with lower numbers of valence electrons. These were generated 
with the program ADPACK, contained in the OpenMX package. We generated fully 
relativistic pseudopotential with the conhguration, he., with 6 valence electrons, 

for Te, and with the 4s^4p^ conhguration for Ge. The pseudopotentials were calculated 
using Troullier-Martins [3l] algorithm with Kleinmann-Bylander factorization [35] . For 
Ge, the cutoff radii for 4s (4p) pseudopotentials were 1.80 (1.70) atomic units (a. u.). 
For Te, the cutoff radii were 2.40 (2.20) a. u. for 5s (5p) pseudopotentials. The radial 
cutoffs for pseudoatomic orbitals were 7.5 and 8.0 a. u. for Ge and Te, respectively. 
Furthermore, for Ge and Te we used the minimal wave functions basis, he., slpl in the 
notation explained in Refs. [36l [37] , which means that in the wave function basis for 
Te and Ge atoms one s and three p orbitals (p^,, Py and p^) were included. The number 
of orbitals was doubled due to the inclusion of electron’s spin. These pseudopotentials 
and the basis set lead to equilibrium lattice parameters at/i=6.0 A, ath = 89.3" and 
Tth = 0.019, which are close to the experimental values. For comparison, the original 
pseudopotentials distributed with OpenMX package give at/i=6.09 A, ath = 89.5° and 
Tth = 0.014. 

For Mn, we used the pseudopotential distributed with OpenMX with 15 valence 
electrons (3s^3p®3(i®4s^ conhguration), together with pseudoatomic basis s?>p2dl. This 
choice increases the number of valence electrons, but since there is only one Mn atom 
in the supercell, the computational time is practically not affected. 

Energy levels of atomic shells with I > 1 are always split by the S-0 interaction. The 
energy difference between split levels, AE = A°i7, where = E{j = 1 + 1/2) — E{j = 
/ — 1/2), is directly related to the S-0 strength. Usually, atomic pseudopotentials for 
j = 1 + 1/2 are calculated assuming that the values of these splittings are those for 
a real atom, and such j-dependent pseudopotentials are called fully relativistic. On 
the other hand, by performing appropriate average [58] over j-dependent pseudopoten¬ 
tials, it is possible to obtain the so-called scalar relativistic pseudopotentials for which 
AE = 0. For these two cases, AE = A^E and AE = 0, the factor fso describing in 
the following the S-0 strength is fso = 1 and fso = 0, respectively. OpenMX package 
offers the possibility to generate j-dependent pseudopotentials for arbitrary fso factors, 
for which AE = fsoA^E. This mimics the variation of the S-0 interaction strength for 
a given shell of an atom. To study the influence of the S-0 interaction on MAE, we 
generated pseudopotentials for Ge and Te for a few values of fso, while for Mn we always 
used the pseudopotential with fso = l, corresponding to the actual atomic S-0 strength. 
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An important feature of OpenMX is the possibility to constrain directions of spin 
polarization for atoms. This is done by adding a properly constructed harmonic potential 
to the total energy functional. To introduce notation which is used below, we briefly 
sketch this construction [39]. In relativistic calculations the wave function is a two 
component spinor 



and the electron density operator h is a 2 x 2 matrix 

n = Y.fu\'4’u><iJu\, ( 2 ) 

where f„ is the occupation number. By projecting this matrix on orbitals of a given 
atom one obtains a matrix N, which can be diagonalized by an unitary matrix U{9, (f) 
dependent on two angles, 6 and (p: 

u[e,^)NuHeM = {^^^ ^J . (3) 

The angles 9 and ip define both the direction of total spin and the difference vr = — 

which is the spin polarization of a given atom. 

To constrain spin of the atom to the direction dehned by angles 9q and 0o, the following 
term is added to the total energy functional: 

E,, = uTr(^(iV-iVo)'), (4) 

where 

= (5) 

If after the selfconsistent procedure the spin direction is the same as the demanded one, 
the energy connected with this term is zero. Indeed, after calculations this energy was 
always negligible compared to other energies of the system, and the final angles 9 and 
if describing the spin directions were practically the same as those demanded provided 
the spin polarization of atoms was not negligible. 

In most calculations, directional constraints were applied only to Mn. Analysis 
of the results shows that the spins of Ge atoms are parallel to the spin of Mn, and 
the directions of Te spins are opposite. Moreover, the calculated MAE does not depend 
signihcantly on whether the constraint was applied to the Mn spin only, or to all atoms in 
the supercell. We conclude therefore that the direction of Mn spin governs the directions 
of Ge and Te spins. 

In all calculations, the Mn spin is rotated in (110) plane and the angle 9 is measured 
in this plane from the [111] direction. Such a rotation corresponds to the situation 
frequently studied in electron ferromagnetic resonance experiments, when the external 
magnetic held rotates in this plane. Gei-j^Mn^jTe on BaF 2 substrates grows along the 
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[111] direction, thus 6 = 0 and 6 = 90° correspond to the Mn spin perpendicular and 
parallel to the plane of growth, i.e., oriented along the [ 111 ] and[ 112 ] crystallographic 
directions, respectively. 

2.3. Other input parameters 

The calculations were performed in part for nonzero densities of free holes p, because 
actual Gei_a;Mn 3 ;Te samples are p-type due to very high concentrations of native defects 
(Ge vacancies). This additional positive charge was compensated by a uniform negative 
background. The hole concentrations assumed in calculations were in the interval 
0 < p < 5 X lO^^cm”^. To speed up the convergence, particularly in the case of nonzero 
hole concentrations, the calculations were performed assuming the electron temperature 
300 K. 

The integrals in numerical calculations are approximated by sums over discrete sets 
of points. Thus, it is important to properly choose the grid for the system in the real 
and the reciprocal space to secure an optimal balance between the accuracy and the 
computational time. Gonvergent results are obtained for the length of the grid in real 
space approximately equal to 0.15 A. Much more important is the number of integration 
points in the Brillouin zone giving convergent results, which, as it is discussed below, 
strongly depends on the concentration of holes. 

2 . 4 . Angular dependence of MAE 

Group theory predicts a general dependence of MAE on the direction of magnetization. 
The knowledge of this dependence enables signihcant saving of computer time. For the 
Oh crystal symmetry, the hrst nonvanishing terms in MAE are of the fourth order in 
magnetization, or in the total spin components Si, i = x,y,z. When the symmetry is 
lowered to for rhombohedral Gei_ 3 ;Mnj.Te, the lowest order terms are proportional 
to S]; and Sf, and MAE can be written as 

Ea{6) = a 2 cos 26 + 62 sin 26 + 04 cos 46 + 64 sin 46 + c, ( 6 ) 

where 6 is the angle between the direction of the Mn spin placed in (1-10) plane and the 
[ 111 ] crystallographic direction, and c is the isotropic part chosen in such a way that 
Ea(6 = 0 ) = 0 . 

In most cases, the rhombohedral angle a = 88° and the displacement r = 0.03 were 
assumed. For these parameters the hrst term in Eq. (jH]) dominates, and MAE may be 
written in the form 

Ba(6) = a2(cos 26-1). (7) 

The remaining terms in Eq. ([S]) are important when the crystal structure is close to the 
perfect cubic, or when two or more Mn ions are placed in the supercell. The validity 
of the htting formula, Eq. ([7]), is fully conhrmed by the obtained numerical results. 
Moreover, the correctness of Eq. ([7]) for rhombohedral Gei_a,Mn 3 ;Te is conhrmed by 
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Figure 1. GeTe band structure for different crystallographic structures: (a) a = 90°, 
T = 0, (b) a = 88°, r = 0, (c) a = 90°, r = 0.03, (d) a = 88°, r = 0.03. The horizontal 
lines at (d) represent positions of the Fermi level for hole concentrations 10^^cm“^, 
3 X 10^^cm“^, and 5 x 10^^cm“^, respectively. 


ferromagnetic resonance measurements [231 ED], which also lead to the conclusion that 
the terms proportional to the fourth order in magnetization are very small. 

3. Factors that determine MAE 
3.1. Band structure of GeTe 

We first discuss the impact of the crystal distortions on the band structure of pure GeTe. 
Figure [D shows the band structure for four sets of structural parameters, namely (a) 
a = 90°, r = 0, (b) a = 88°, r = 0, (c) a = 90°, r = 0.03, and (d) a = 88°, r = 0.03. 
In the first case the crystal symmetry is Oh, and it is lowered to for (b), and to 
for both (c) and (d). 

When the symmetry is lower than Oh, the T and L points in the Brillouin zone are 
nonequivalent. They are defined by T = |(bi + b 2 + bs) and L = ^bi, respectively, 
where bi, b 2 , bs are the primitive vectors of the reciprocal lattice. Other points in the 
Brillouin zone of a rhombohedral lattice are defined in Ref. [TT] , 

As it follows from Fig. (H the impact of the rhombohedral distortion on the band 
structure is much weaker than that of the internal displacement of the Ge and Te 
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sublattices, since the latter opens the band gap of about 0.5 eV at both L and T points. 
One can also observe that the decrease of the angle a leads to the decrease of band 
energies in the vicinity of the L point compared to those in the vicinity of T, while the 
increase of parameter r has the opposite effect. 

The band structures shown in Fig. 1 are calculated for wavefunction basis slpl. 
Qualitatively, the same results are obtained for enlarged bases s2p2 and s3p3. 

It is known that the LDA calculations underestimate the energy gap. Due to this 
effect, GeTe in the Oh structure is found to be semimetallic with zero energy gap. This 
is also the case when the rhombohedral distortion is accounted for but the internal 
displacement r is neglected. Fig. [H For the present work it is important that the 
system has a non-vanishing energy gap for the used crystal parameters, close to the 
experimental ones, a = 88° and r = 0.03. This enables a comparison of magnetic 
anisotropy in insulating and conducting crystals, and an analysis of the role of free 
carriers. For hole concentrations of the order of 10^^ cm“^, typical for Gei_ 3 ;Mn 2 ,Te 
[23], the Fermi level lies deep in the valence band and the precise band structure in the 
vicinity of the band maximum is not important. 

3.2. Dependence of MAE on the hole concentration 

Figure E] presents the angular dependence of MAE for the 2x2x2 supercell containing 
one Mn ion for 3 hole concentrations. The calculated values are shown as points, and 
the curves are obtained by htting the results according to Eq. ((Tj), i.e., using one htting 
parameter 02 . In all cases, the angular dependencies of are symmetric relative to 
6 = 0 . 

The values of the coefficient 02 , calculated for several hole concentrations p ranging 
from zero to 5 x 10^^ cm“^, are summarized in Fig. [31 The obtained results show that 
both the absolute value of MAE and its sign depend on the concentration of free carriers. 
Indeed, the presence of free holes can enhance the magnitude of MAE by more than 
one order of magnitude compared to the insulating case. As discussed below, this is 
caused by a larger and more spatially extended spin polarization around Mn in metallic 
samples. 

Regarding the sign of MAE, three ranges of hole concentrations can be 
distinguished. First, for vanishing or very small p, the coefficient 02 is positive, i.e., 
the equilibrium orientation of Mn spin is perpendicular to the [111] direction. Second, 
in the window 1 x 10^® cm“^ - 2.5 x 10^^ cm“^, 02 is negative. Finally, for higher p, 
02 is positive again. Moreover, with the increasing hole concentration a non-monotonic 
behavior is observed. 

The impact of free carriers on magnetic properties reflects pecularities of the 
band structure. Ghanges of the direction of easy axis caused by changes in carrier 
concentration were observed and analyzed for Gai-ajMn^-As [T31 jUj. Similarly, 
in Sni_a;Mna;Te, the cubic anisotropy constant depends on the hole concentration 
|2Uj . Another effect related to the changes in carrier concentration was observed in 
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Pbi_j^_ySna,MnyTe diluted magnetic semiconductor, in which the Curie temperature 
depends non-monotonically on the hole concentration |33] . The effect was explained by 
signihcant modihcations of the RKKY interaction caused by the two-band structure {L 
and S bands) of the valence band |15]. The multiband structure of the valence band 
is characteristic for IV-VI compounds, and it is absent in e.g. Gai_j;Mnj;As. In the 
case under consideration, the number of valence bands occupied with holes depends 
on the Fermi energy, Ep. In Fig. [T]i we show the Fermi level position for three hole 
concentrations. For the highest concentrations, as much as four bands are occupied 
(at T, L, and those on both X-T and T-U directions). Each of those bands has an 
its own characteristic orbital composition, strength of the S-0 coupling, ^f-factor, and 
coupling with Mn. For this reason, the non-monotonic dependence of MAE on the hole 
concentrations, shown in Fig. [3l can be related to progressive occupation of consecutive 
valence bands. 

Figure m presents the results of calculations for the 3x3x3 supercell containing 
one Mn ion for both the insulating and the metallic crystal with p = 5 x 10^^ cm“^. As 
it follows from the hgure, the directions of easy axes are the same as for the 2x2x2 
supercell. However, the magnitude of 02 is about twice higher for the 3x3x3 supercell 
in both insulating and metallic cases. These results show that the Mn ions in the 2x2x2 
supercell cannot be treated as isolated, even in the insulating crystal. 

Finally, also shown in Figs. ElandlUis the convergence with respect to the number 
of integration points (NIPs) in the Brillouin zone. The convergence is easier to achieve 
for the insulating system than for the metallic one. In the former case, the obtained 
dependencies are practically the same for three consecutive NIPs, 7^, 8^ and 9^. When 
the hole concentration is hnite, the convergence is much slower. Unfortunately, due to 
limitations in computational resources we could not perform calculations for NIP greater 
than 14^ for the 2x2x2 supercell for all the considered cases. However, for one selected 
concentration p = 2.5 x 10^° cm“^, NIP was increased to 16^=4096 and 18^=5832, and 
the results shown in Fig. 2b confirm that a satisfactory accuracy of about 5 per cent 
was achieved. 

The increase of the supercell size leads to the so-called folding of the Brillouin zone. 
Due to the folding, the grid spacing obtained with NIP=12^ used for the 2x2x2 su¬ 
percell is the same as NIP=8^ used for the 3x3x3 supercell, and both correspond to 
the grid generated with NIP=24^ for the simple cubic unit cell. From Figs [2] and 0] it 
follows that the accuracy of the coefficient 02 calculated for metallic systems is about 
ten per cent, and less than one per cent for insulating systems. 


3.3. Dependence of MAE on the macroscopic structural parameters a and r 

Figure [5] shows the angular dependence of MAE on crystal structure parameters a and 
r for two concentrations, p = 0 and p = 10^^ cm“^. The impact of the lattice parameters 
is qualitatively different in the insulating and conducting crystal. The change in the 
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Angle (deg) 

Figure 2. Angular dependence of energy of magnetic anisotropy for the 2x2x2 
supercell and three different hole concentrations: (a) p=0, (b) p = 2.5 x 10^° cm“^, 
and (c) p = 5 X 10^^ cm“^. The numbers in legends are the numbers of grid points in 
the Brillouin zone. 

internal displacement from r = 0.03 to r = 0.02 (continuous and dotted lines in Fig. 0 
respectively) results in significant change of the amplitude of MAE in the insulating 
case, while in the conducting case it is less important. On the other hand, changes 
in the rhombohedral angle a, which are not very important in the insulating crystal, 
produce much bigger effect in the conducting case. Similar differences are discussed in 
the next subsection, they are related to the fact that the spatial region around the Mn 
ion contributing to MAE is more localized in the insulating case. 

Let us also remind that Gei_a;Mn 3 ;Te containing more than 20 at. % of manganese 
is not rhombohedral but cubic, and the easy axis of magnetization lies in plane of the 
layer (i.e., it is perpendicular to the [111] direction). Although our calculations were 
performed for much lower manganese concentration and quantitative conclusions cannot 
be drawn, the direction of changes in MAE values is in accordance with experiment. 

3.4- Dependence of MAE on the geometry of the Mn surrounding 

In the metallic cases, the calculated amplitude of MAE is at least an order of magnitude 
larger than that observed in experiment. However, one should keep in mind that the 
calculations were done for supercells which represent a perfect translationally invariant 
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Figure 3. The dependence of the coefficient 02 of MAE (Eq. 2) on the hole 
concentration. The solid line is to guide the eye only. For p > 0 all the calculations 
were performed for 2744 integration points in Brillouin zone. 


infinite lattice. As it was mentioned in the Introduction, in real mixed crystals there 
is always a microscopic disorder, as a result of which each Mn ion is in a different 
surrounding. Moreover, the experimental Mn concentrations are an order of magnitude 
higher than those considered here, and thus the Mn-Mn coupling must be taken into 
account. 

Figure [H] shows the angular dependence of MAE for 3 atomic conhgurations of the 
Mn neighbours. In all cases, the symmetry of the crystal is rhombohedral. The two 
left panels. Figs. Et'' andEb, show MAE for the insulating case, and Fig. Et for the hole 
concentration p = 5 x 10^^ cm“^. 

First, the continuous lines represent MAE for the ideal lattices, characterized by 
a = 88° and r=0.03, and they already were presented in Fig. [2J Due to the non¬ 
vanishing internal displacement r, 3 out of 6 Mn-Te distances are equal 3.189 A, and 
the 3 remaining are 2.819 A. 

Second, the dotted lines show the results for the conhguration when the Mn ion is 
moved along the main diagonal of the supercell by ao\/3r, and all distances to the six 
nearest neighbouring Te ions are the same. In the metallic case, the displacement of the 
Mn ion has a very small influence on MAE, while in the insulating case the impact is 
substantial. More specifically, the absolute values of changes of Ea{6 = —90°) are 0.2 K 
and 0.8 K in the metallic and insulating cases, respectively, while the relative changes 
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Angle (deg) 


Angle (deg) 


Figure 4. Angular dependence of energy of magnetic anisotropy for the 3x3x3 
supercell for two different hole concentrations: (a) p=0 and (b) p = 5 x 10^^ cm“^. 
The numbers in legends are the numbers of grid points in the Brillouin zone. 


are 3% and 277%, respectively. 

This difference is connected with the much larger spatial range and magnitude of 
spin polarization of valence electrons in the metallic than in the insulating case. This 
issue is discussed in detail in the following. Consequently, in the metallic sample the 
spin polarization extends over the whole supercell, and the very localized perturbation 
of the lattice caused by a displacement of a single ion has a small influence on MAE. 
On the other hand, in the insulating case the range of spin polarization is much smaller, 
and therefore the system is sensitive to the perturbation in the nearest neighbourhood 
of Mn. That is why the changes of MAE are also larger. 

Finally, the broken lines display the results obtained for conhgurations with 
partially optimized atomic positions. In the optimization procedure, the Mn ion together 
with its nearest 12 Ge and 6 Te neighbours were allowed to move, while the remaining 
atoms in the supercell were kept at their initial ideal positions. The displacements of 19 
atoms lead to a much larger perturbation of the lattice than the displacement of the Mn 
atom only, and the changes of MAE are signihcant in both metallic and insulating cases. 
For this partially relaxed Gei_ 2 ;Mna;Te lattice, the Mn-Te bond lengths are 3.043 A and 
2.796 A. Thus, although the replacement of Ge by Mn atom does not lead to such a large 
local lattice deformation as in other IV-VI semiconductors (for example, in Pbi-ajMuj-Te 
the difference between ionic radii of Pb and Mn is much larger than that between Ge 
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Figure 5. Angular dependence of MAE for the 2x2x2 cell for (a) p = 0, NIP=512 
and (b) for p = lO^^cm"^, NIP=2744 for different crystal structure parameters a and 

T. 

and Mn), the impact of local distorsions on MAE is certainly not negligible. 

The details of MAE for p = 0 and the partially relaxed lattice are shown in the 
lower left panel of Fig. El The broken line is a £t to the calcnlated points according to 
formula resulting from group theory considerations 

Ea = a 2 cos 26 + a^c cos^ 6 {cos^O — 7/Qj + 04 ^ sin^ 6 cos 6 + c, ( 8 ) 

where the second and the third term on the right hand side are proportional to the an¬ 
gular dependent parts of the Legendre polynomials P 4 and Pf, respectively. Only such 
terms are allowed for rhombohedral symmetry if we consider terms up to fourth order 
in spin components. One can observe that such a fitting reproduces well the calculated 
values. 




3.5. MAE of the Mn-Mn pairs 

We also performed preliminary calculations for 2x2x2 supercells containing two Mn 
atoms. For such supercells there are 21 possible nonequivalent conhgurations. For most 
of them the symmetry is not rhombohedral any more, thus Eq. ([7]) is no longer valid, 
and the term 62 sin 20 must be added to the fitting formula for Ea. The coefficients 02 
and 62 depend on both the configuration and the hole concentration. Interestingly, the 







Calculations of magnetic anisotropy energy of Gei-xMuxTe 


15 



Figure 6. Impact of the lattice parameters on the angular dependence of MAE for 
the 2x2x2 cell for (a), (b) p = 0, NIP=512 and for (c) p = 5 x 10^^cm“^, NIP=2744. 
The meaning of continuous, dotted and broken lines is explained in the text. 


average of 62 over configurations almost disappears, indicating that although for most 
conhgurations the crystal symmetry is not rhombohedral, after the averaging procedure 
the functional dependence of MAE is in accordance with the macroscopic rhombohedral 
symmetry. 

The situation is similar to that considered in Ref. [I 6 ], where the single ion 
anisotropy of Mn in Sni_a;Mn 3 .Te crystal was analyzed. In this work, the chemical 
microscopic disorder caused by the random placement of Mn ions in the lattice was 
taken into account. For each Mn ion, the amplitude of MAE was much larger than 
that observed in experiment, and the angular dependence did not correspond to the 
Oh symmetry because the nearest neighbourhood of Mn ions had not such symmetry. 
However, after the averaging procedure the average MAE was reasonably close to the 
experimental values, and the angular dependence was consistent with the macroscopic 
cubic symmetry of the crystal. 
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4. Mechanism of magnetic anisotropy 
4 . 1 . Spin polarization 

Atomic spin polarizations tt, defined in Section II, was analyzed for the 3x3x3 supercells 
for two hole concentrations, p = 0 and p = 5 x 10^^ cm“^. The polarization depends on 
the atomic species, the distance from Mn, and on the hole concentration. 

For p = 0, the Mn polarization = 4.76, which is close to = 5 for the free 
atom. The values of and are much smaller, of the order of 0.01, and practically 
vanish (i.e., they are less than 0.001) beyond the first Te and Ge neighbours. 

In the metallic case, the spin polarization of Mn depends on the hole concentration, 
and changes from 4.76 to 4.68 for the considered concentration range. The main effect 
induced by the presence of free holes is the pronounced increase of of the nearest Mn 
neighbours from 0.01 to 0.07. Also the spatial range of spin polarization of Te anions 
is larger, since it vanishes only for distances greater than 10 A from Mn. On the other 
hand, spin polarization of Ge cations is similar to that in the insulating case. These 
results are explained by observing that in the metallic phase, Mn polarizes free holes, 
i.e., the states from the valence band top. These states are mainly built up from the 
atomic states of Te, which are polarized by Mn. The contribution of Ge orbitals to the 
hole states is much smaller, and so is the spin polarization of Ge. As it is shown in 
the next subsections, the strength of the spin polarization and its range are the decisive 
factors that determine MAE. 

4-2. Dependence on the spin-orbit strength 

Figure [7] displays the angular dependence of MAE for different values of fso, the factor 
that tunes the spin-orbit interaction strength for Ge and Te. The figure compares the 
results for the insulating p = 0 case. Fig. [Tti, with those for the metallic one with 
p = 0.5 X lO^^cm”^, Fig. [Tb- 

When fso = 0, i. e., the S-0 coupling is non-vanishing for the Mn ion only, the 
angular dependencies of MAE are nearly the same in the insulating and the metallic 
cases (note that the energy scales for left and right panels of Fig. [7|are different), and 
the anisotropy is very weak. It follows from the non-vanishing S-0 coupling of Mn, 
together with the hybridization between its orbitals and the orbitals of neighbouring 
atoms. 

When the S-0 interaction is turned on, the neighbouring atoms contribute to MAE 
as well. The reason is that for a given crystal structure, from the total energy point 
of view, each atom has its preferred direction of the polarization which, in general, is 
not the same as the actual direction of the manganese spin. That is why a change of 
direction of Mn spin changes total energy. 

The impact of the increasing strength of the S-0 interaction is different in the 
insulating and the conducting systems. In the former case, the anisotropy changes sign, 
but MAE remains small. In the metallic regime, the amplitude of MAE can increase 
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Figure 7. Angular dependence of energy of magnetic anisotropy for a 2 x 2 x 2 supercell 
with a = 88°, t = 0.03 containing 1 Mn ion for hole concentration (a) p = 0 and (b) 
p = 0.5 X 10^^cm“^. The respective curves correspond to different S-0 strength for Ge 
and Te pseudopotentials. 


by an order of magnitude or more, and it can also change sign depending on the actual 
hole concentration, see Fig. 3. We also found that the polarizations of Ge and Te 
atoms and their directions do not depend on fso signihcantly. Actually, the differences 
in polarizations for Ge or Te for different fso are less than 0.005. 

These results clearly show that even in the insulating case the magnetic anisotropy 
is not related to Mn ion only but also to atoms in the Mn neighbourhood, the spatial 
range of which depends on the hole concentration. 


4-3. Components of the total energy 

Total crystal energy can be decomposed into the contributions of kinetic energy, 
electrostatic electron-electron and electron-core interactions and exchange-correlation 
energy. In Fig. Owe show changes of the total energy, changes of its components, 
Ekim Eel, Exc, and also changes of the band energy, Eband, i- e. the sum of Kohn - 
Sham hamiltonian eigenvalues over occupied states. We note that the electrostatic term 
contains not only electron - electron and electron - core interactions, but also the S-0 
interaction, because the pseudopotentials are spin-dependent. From Fig. [H] it follows 
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Figure 8. Anglular dependence of energy of magnetic anisotropy components for a 
2x2x2 supercell, the hole concentration p = 2.5 x 10^^cm“^, a = 88°, and r = 0.03. 


that the changes of components of the total energy are much larger than the changes of 
which is typical for DFT calculations. The results of all the considered cases are 
qualitatively similar to those shown in Fig. [HI They clearly indicate that the changes 
of spin direction of the Mn ion are accompanied by spatial redistribution of electron 
density. 

f.f. Microscopic picture of magnetic anisotropy 

The above considerations allow us to propose the following picture for the mechanism of 
magnetic anisotropy. Assume that we know the total energies of the system, Ei and E 2 , 
for two directions of the Mn spin, hi and h 2 , respectively. As it was discussed above, 
the direction of the Mn spin determines the spin directions of remaining atoms in the 
system. 

Let us assume now that we know the solutions to the Kohn-Sham problem for the hrst 
direction of the Mn spin, hi, and let us denote this set of eigenfunctions, corresponding 
to El, by {'01 }• Let us now change, using a unitary transformation, the directions of 
spins in {'0i} in such a way that after the rotation the polarization of atoms correspond 
to the second direction of Mn spin, h 2 . The total energy E[, calculated with the new 
set of wave functions {"01}, in general, is different from Ei. This unitary transformation 
does not change electron densities and (which are the eigenvalues of electron 
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density operator at a given point, Eq. ([2])) thns it does not change the kinetic, electron- 
electron, electrostatic part of the electron-core, and the exchange correlation energies. 
The only term in the Kohn-Sham hamiltonian which does not commnte with the 
transformation is that describing the electron-core interaction, becanse pseudopotentials 
are spin dependent. Thus, the difference between Ei and E[ is solely due to spin-orbit 
interaction. 

If we now take this new set of the wave functions, as a starting point for the 

selfconsistent procedure, at the end we obtain total energy E 2 and the set of wave 
functions {'^> 2 }- The directions of polarization of atoms calculated for {' 1 /^ 2 } are the 
same as for however the resulting spatial electron density change. 

Thus the change of spin directions is accompanied by changes in spatial electron 
distribution. These changes are, of course, very small and may be neglected in 
phenomena in which we are interested in energies of the order of 1 eV, however they are 
critical for the present problem of magnetic anisotropy. One should also notice that the 
spatial redistribution of electrons indicates that the so-called magnetic force theorem 
m, sometimes used in calculations of magnetic anisotropy [ISl SS] is not applicable 
in the present case. Selfconsistent calculations are necessary for each direction of the 
manganese spin. 

On the basis of the presented results we can state that the anisotropy in 
Gei_a:Mna;Te, particularly in the conducting case, is certainly not a single Mn ion 
anisotropy. This is a collective effect in the sense that all the atoms in the supercell 
give contribution to the phenomenon. For higher manganese content, as is in the cases 
studied in experiment, this is also the collective effect. This conclusion is in accordance 
with the Zener model [a [3] where the magnetic anisotropy is also due to free carriers. 
However, ab initio calculations show that the differences in occupancies of band levels 
are not the only cause, the redistribution of electron densities should be also taken into 
account. 

5. Summary 

Magnetic anisotropy energy MAE in the diluted ferromagnetic semiconductor 
Gei_a,Mna;Te was calculated within the Local Spin Density Approximation. This 
approximation is shown to be accurate enough to describe MAE effects of the order 
of 0.1 Kelvin. Low concentrations of Mn were considered, in which magnetic ions are 
interacting only weakly. The analysis included the impact of factors that determine 
MAE, namely (i) the concentration of free holes, (ii) the microscopic conhguration of 
atoms in the surrounding of Mn, (iii) the macroscopic symmetry of the crystal (which 
in the case of Gei_ 2 ;Mna,Te can be either cubic or rhombohedral), and (iv) the strength 
of the S-0 interaction in valence band. The main conclusions are as follows. 

First, magnetic anisotropy is largely determined by the presence of valence band 
holes. The increase of the hole concentration leads to both the variation of the MAE 
strength and to the change of its sign, i.e., to the reorientation of the magnetization 
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easy axis from parallel to perpendicular to the [111] direction. In particular, for hole 
concentrations of the order lO^^cm”^, typical of Gei_a;Mna.Te, the calculated direction 
of the easy axis of magnetization agrees with that observed in recent experiments 
on Gei_j,Mn 3 ;Te epitaxial layers deposited on BaF 2 (111) substrates (perpendicular 
magnetic anisotropy) [23l EHl EHl EQ] . 

Second, MAE is surprizingly sensitive to both the macroscopic crystal structure and 
the atomic conhguration of the Mn neighbourhood. For example, a change of the crystal 
parameters from a = 88° and r = 0.03 to a = 89° and r = 0.02 reduces MAE by as 
much as one order of magnitude. A comparable sensitivity is found for the dependence 
of MAE on the local atomic conhguration. 

Third, the calculated angular dependencies of MAE fully agree with those predicted 
by group theory. For Gei_a;Mn 3 ;Te in the rhombohedral phase, the terms of the 
second order in Mn spin components satisfactorily describe the calculated results. This 
dominance of the quadratic terms is consistent with the experimental results [231 ITU] . 
On the other hand, higher order terms cannot be neglected for cubic Gei_a;Mn 3 ;Te. 

Finally, the role of the spin-orbit coupling was elucidated by calculating MAE as a 
function of the S-0 strength of the band states of host Gei_a;Mna;Te crystal. The S-0 
of the host is particularly important in conducting samples. 

The obtained results are related to the fact that the spin polarization of the host 
atoms have a decisive impact on MAE. Since the states from the top of the valence 
bands are mainly composed from the p orbitals of Te, they determine spin polarization 
of the valence states. The insulating and the metallic case differ by the spatial extension 
of the spin polarization: in the former case only the hrst two coordination shells are spin 
polarized, while in the latter case spin polarization extends over several lattice constants. 
This varying degree of localization explains a number of the calculated trends in MAE. 
Moreover, in the presence of S-0 interaction the changes in the Mn spin direction induce 
changes in the electron density, the impact of which on the total energy was discussed 
in the last section. 

Regarding the comparison with experiment one notices that a description of MAE 
in real samples containing 10-20 per cent of Mn must include Mn-Mn magnetic coupling, 
much stronger in metallic than in insulating samples. Moreover, incorporation of the 
microscopic chemical disorder, which implies in particular local deformations of the 
lattice in the vicinity of Mn, is necessary. On the other hand, the orientation of 
easy axis (perpendicular magnetic anisotropy) observed in Gei_a:Mna;Te/BaF 2 (HI) 
layer is already explained by the present results. Although the calculations were 
performed for Gei_ 3 ;Mn 3 ;Te, they directly apply to closely related IV-VI ferromagnetic 
semiconductors, Sni_a,Mn 3 ;Te and Pbi_ 3 ;_ySn 3 ;Mnj^Te. The proposed microscopic 
mechanism of magnetic anisotropy should also be applicable to other groups of diluted 
ferromagnetic semiconductors. 
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